Journal  of  Power  Sources  262  (2014)  147-161 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


Calibration  method  for  carbon  dioxide  sensors  to  investigate  direct 
methanol  fuel  cell  efficiency 

M.  Stahler*,  A.  Burdzik 

Institute  of  Energy  and  Climate  Research,  (IEK-3:  Electrochemical  Process  Engineering),  Forschungszentrum  Jiilich  GmbH,  Wilhelm- Johnen-StraJSe, 
52425  Jiilich,  Germany 


CrossMark 


HIGHLIGHTS 
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•  The  accuracy  is  improved  by  a  factor  of  2—30  —  depending  on  the  metering  range. 
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Methanol  crossover  is  a  process  in  direct  methanol  fuel  cells  which  causes  significant  reduction  of  cell 
efficiency.  Methanol  permeates  through  the  membrane  electrode  assembly  and  reacts  at  the  cathode 
with  oxygen  to  form  carbon  dioxide.  This  process  is  undesirable  because  it  does  not  generate  electric 
energy,  but  rather  only  increases  heat  production. 

Different  procedures  have  been  used  for  the  investigation  of  this  crossover.  One  method  uses  the 
detection  of  carbon  dioxide  in  the  exhaust  gas  of  the  cathode  by  means  of  a  carbon  dioxide  sensor. 
This  technique  is  inexpensive  and  enables  real-time  measurements  but  its  disadvantage  is  the  low 
accuracy. 

This  paper  demonstrates  a  simple  method  to  generate  gas  mixtures  for  the  calibration  of  the  sensor  in 
order  to  increase  the  accuracy.  The  advantages  of  this  technique  consist  in  the  fact  that  only  the  existing 
devices  of  a  direct  methanol  fuel  cell  test  rig  are  needed  and  that  the  operator  can  adjust  the  carbon 
dioxide  concentration  for  the  calibration  process.  This  is  important  for  dealing  with  nonlinearities  of  the 
sensor.  A  detailed  error  analysis  accompanies  the  experiments.  At  the  end  it  is  shown  that  the  accuracy  of 
the  determined  Faraday  efficiency  can  be  improved  by  using  the  presented  calibration  technique. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Direct  methanol  fuel  cells  (DMFC)  are  interesting  for  supplying 
electrical  energy  in  applications  with  power  requirements  of  up  to 
a  few  kilowatts.  One  advantage  of  a  DMFC  system  in  comparison  to 
other  fuel  cell  systems  is  the  high  energy  density  of  the  methanol 
which  enables  a  long  operating  time.  A  disadvantage  of  the  DMFC  is 
the  methanol  permeation  from  the  anode  of  a  membrane  electrode 
assembly  (MEA)  to  the  cathode.  On  the  cathode  side,  the  methanol 
reacts  with  oxygen  to  form  carbon  dioxide.  This  crossover  reduces 
the  fuel  utilization  and  therefore  the  efficiency  of  the  system. 
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Researchers  employ  different  techniques  to  investigate  the 
correlations  between  operating  conditions,  materials  used  for  the 
manufacturing  of  the  MEA,  and  methanol  permeation.  Some  of 
them  analyze  the  cell  voltage  [1,2],  whereas  others  investigate  the 
exhaust  gas  composition  of  the  DMFC  [3-7]  to  determine  the 
strength  of  the  crossover.  While  precise  techniques  for  in¬ 
vestigations  of  the  gas  composition,  such  as  mass  spectroscopy  and 
gas  chromatography,  are  expensive  and  time-consuming,  a  cheaper 
infrared  carbon  dioxide  sensor  also  enables  real-time  measure¬ 
ments.  Unfortunately,  the  accuracy  of  this  type  of  sensor  depends 
strongly  on  the  metering  range.  As  a  result,  the  relative  uncertainty 
of  the  measurement  can  vary  between  2%  and  more  than  20% 
within  the  measuring  range. 

For  DMFC  investigations,  operating  conditions  must  be  varied  in 
order  to  analyze  the  correlations.  During  these  variations,  the  car¬ 
bon  dioxide  concentration  in  the  exhaust  gas  can  differ  very 
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strongly,  depending  on  the  operation  mode.  Several  approaches 
can  be  employed  to  avoid  strong  differences  in  the  uncertainties  of 
these  measurements.  First  of  all,  it  is  possible  to  use  several  sensors 
with  different  metering  ranges  simultaneously.  This  circuitous 
method  reduces  the  advantage  of  a  cost-efficient  measurement 
device.  One  further  option  is  to  use  commercial  calibration  gases  for 
a  precise  calibration  of  a  sensor  within  a  wide  metering  range.  The 
problem  with  this  method  is  that  a  nonlinearity  in  the  metering 
range  of  the  sensor  requires  many  calibration  gases  with  different 
concentrations  to  achieve  high  calibration  accuracy.  However, 
before  calibration,  the  nonlinear  metering  ranges  of  a  sensor  are 
unknown  so  that  the  right  choice  of  the  necessary  calibration  gases 
is  difficult. 

The  work  presented  here  focuses  on  the  development  of  an 
easily  performable  technique  to  calibrate  a  carbon  dioxide  sensor  in 
a  DMFC  test  rig  within  a  wide  metering  range.  The  relation  to  the 
experimental  station  means  that  all  relevant  components  for  the 
calibration  process  can  be  taken  from  a  DMFC  test  rig.  In  this  study, 
the  desired  relative  accuracy  for  a  metering  range  of  0.2%  v/v  to 
20.0%  v/v  is  less  than  3%.  The  idea  of  the  method  is  to  generate 
carbon  dioxide  by  an  electrochemical  reaction.  If  the  reaction  is  free 
of  side  reactions,  the  strength  of  the  generation  can  be  determined 
accurately  by  measuring  the  electric  current  between  the  elec¬ 
trodes.  The  continuously  produced  carbon  dioxide  is  mixed  with  a 
nitrogen  flow,  whereat  the  flow  rate  of  the  nitrogen  can  be 
measured  precisely  by  a  laminar  flow  element.  This  method  allows 
the  production  of  calibration  gas  with  an  adjustable  carbon  dioxide 
concentration.  In  this  way,  the  nonlinear  metering  ranges  of  the 
sensor  can  be  analyzed  very  accurately  and  a  wide  range  sensor 
calibration  is  possible. 

By  means  of  the  improved  sensor  calibration,  the  accuracy  of  the 
Faraday  efficiency  from  direct  methanol  fuel  cell  investigations  can 
be  determined  more  precisely.  A  reduction  of  the  uncertainty  of  the 
calculated  Faraday  efficiency  by  a  factor  of  up  to  4  is  demonstrated. 

2.  Experimental 

The  following  sections  present  the  experimental  setup  used  for 
the  investigation  of  DMFC  MEAs.  After  a  short  description  of  the 
devices  four  different  operating  modes  are  explained.  One  mode  is 
used  to  calibrate  the  carbon  dioxide  sensors,  the  other  modes  are 


important  for  DMFC  operation.  The  accuracies  of  the  devices  are 
mentioned  and  form  the  basis  for  the  following  error  analysis. 
Finally,  the  fault  effect  analysis  is  complemented  by  some  necessary 
basic  inspections. 

2.1.  Experimental  setup 

In  the  middle  of  Fig.  1,  which  is  divided  in  side  A  and  side  B,  an 
MEA  is  located  between  two  flow  fields  and  two  end  plates.  On  the 
left  a  peristaltic  pump  (Reglo  digital,  Ismatec)  is  used  to  pump 
aqueous  solution  of  methanol  or  formic  acid  from  a  tank  to  a  pre¬ 
heater  which  tempers  the  solution  to  operating  temperature.  Af¬ 
terward,  a  gas  separator  removes  the  dissolved  gas  and  the  solution 
flows  into  the  test  cell.  During  DMFC  operation  methanol  is 
oxidized  at  the  anode  (side  A)  forming  carbon  dioxide  which  is 
mixed  with  the  aqueous  methanol  solution. 

At  the  anode  output  a  gas  separator  is  used  to  remove  carbon 
dioxide  from  the  outlet  solution.  The  carbon  dioxide  is  separated 
from  the  liquid  phase  by  injecting  nitrogen  gas  bubbles  into  the 
liquid  phase  via  a  porous  glass  plate.  The  nitrogen  volume  flow  is 
controlled  by  a  mass  flow  controller  (5850S,  BROOKS,  adjustment 
range:  0-1  slpm  (standard  liter  per  minute),  accuracy:  0.4% • 
metering  range  +  0.2%  reading,  all  flow  rates  in  this  study  are 
related  to  standard  conditions:  ps  =  101,325  Pa,  Ts  =  273.15  K).  At 
the  end,  the  liquid  phase  flows  through  a  valve  back  to  the  tank  or 
to  the  waste  water  system,  depending  on  the  operation  mode. 

After  the  separator  the  water  vapor  of  the  gas  output  is 
condensed  out  at  room  temperature.  The  gas  temperature  is 
recorded  by  a  digital  voltmeter  (Hexagon  720,  BEHA,  temperature 
accuracy:  0.4  °C)  and  the  absolute  pressure  is  determined  by  a 
pressure  transmitter  (PTU30T,  Vaisala,  pressure  accuracy:  45  Pa). 
The  carbon  dioxide  volume  concentration  crA  in  the  gas  mixture  is 
quantified  by  a  carbon  dioxide  sensor  (GMT221  -  transmitter, 
GMP221  —  sensor,  Vaisala,  metering  range:  20%  v/v). 

At  side  B  of  Fig.  1  a  nitrogen  or  air  gas  flow,  controlled  by  a  mass 
flow  controller  (5850S,  BROOKS,  adjustment  ranges:  0-1  slpm  or 
0-0.12  slpm,  accuracy:  0.4% •  metering  range  +  0.2%  reading)  can 
be  humidified  and  preheated  before  entering  the  test  cell.  To 
measure  the  relative  humidity,  the  output  gas  at  side  B  can  be 
preheated  up  to  130  °C.  This  temperature  is  chosen  because  the 
relative  humidity  of  the  gas  is  smaller  and  therefore  the  sensitivity 
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Fig.  1.  Illustration  of  the  experimental  setup  used  for  the  investigations. 
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of  the  humidity  sensor  is  higher  in  comparison  to  the  operating 
temperature  of  the  test  cell.  After  the  humidity  measurement,  the 
water  of  the  gas  is  condensed  out  at  room  temperature.  In  order  to 
determine  the  temperature  TB,  the  pressure  pB,  and  the  carbon  di¬ 
oxide  volume  concentration  crB,  the  same  device  types  are  used  as 
described  for  side  A.  The  metering  range  of  the  carbon  dioxide 
sensor  is  5%  v/v. 

To  adjust  the  electric  current  between  side  A  and  side  B,  an 
electronic  load  (ZS512,  H&H)  is  used  in  combination  with  a  power 
adapter  (ES015,  Delta  Elektronika).  The  electric  current  I  is 
measured  by  a  multimeter  (Hexagon  720,  Beha,  accuracy:  1%- 
metering  range  +  0.1%  reading).  For  electric  currents  smaller  than 
0.5  A  a  metering  range  of  0.5  A  is  used.  For  larger  currents  a 
metering  range  of  10  A  is  employed. 

The  test  cell  consists  of  single  meander-shaped  anode  and 
cathode  flow  fields  in  combination  with  an  MEA  produced  in  house. 
The  catalyst  loading  at  side  A  is  2  mg  Pt/Ru  per  cm2  (HISPEC  12100, 
Pt/RU/C,  Johnson  Matthey),  the  catalyst  loading  at  side  B  is 
2  mg  Pt  per  cm2  (HISPEC  9100,  Pt/C,  Johnson  Matthey),  and  the 
quadratic  cell  area  is  17.6  cm2.  The  thickness  of  the  MEA  is  600  pm. 
A  membrane  of  type  Nation®  115  is  used  for  the  MEA  preparation. 

2.2.  Operating  modes 

The  setup  can  be  operated  in  different  modes  which  are 
described  in  this  section.  Depending  on  the  research  question  four 
different  modes  are  important:  the  DMFC  mode,  the  pseudo-half¬ 
cell  mode  (phc  mode),  the  reverse  mode,  and  the  calibration  mode. 

2.2.1.  DMFC  operating  mode 

In  this  operating  mode  the  tank  in  Fig.  1  is  filled  with  aqueous 
methanol  solution.  The  solution  is  pumped  into  side  A  of  the  cell 
(anode)  where  different  processes  run  in  parallel,  see  Fig.  2X. 

If  an  electric  current  is  adjusted  by  the  power  source,  the 
methanol  is  oxidized  at  the  anode.  The  generated  protons  flow 
through  the  membrane  to  side  B  (cathode),  the  electrons  get  to  the 
cathode  via  the  external  power  source,  and  together  with  oxygen 
water  is  formed.  This  is  process  1  in  Fig.  2X. 

Irrespective  of  whether  an  external  current  is  adjusted  or  not, 
the  methanol  permeates  through  the  MEA  from  anode  to  cathode. 
The  oxygen  at  the  cathode  side  oxidizes  the  methanol  forming 


carbon  dioxide  and  water.  This  undesired  methanol  permeation 
process  does  not  generate  electric  energy,  but  rather  only  increases 
heat  production.  This  is  process  3  in  Fig.  2X. 

A  common  approach  for  investigating  the  methanol  permeation 
is  to  measure  the  carbon  dioxide  in  the  exhaust  gas  of  the  cathode. 
For  this  approach,  it  is  important  to  know  how  strong  carbon  di¬ 
oxide  diffuses  from  the  anode  to  the  cathode  because  in  this  case 
(process  2  in  Fig.  2X)  the  carbon  dioxide  in  the  exhaust  air  of  the 
cathode  does  not  correlate  with  the  methanol  permeation  process. 

In  case  of  DMFC  operating  mode  the  output  at  the  anode  con¬ 
sists  of  water,  carbon  dioxide,  and  methanol.  The  valve  at  the  anode 
directs  the  liquid  phase  to  the  waste  water  system,  the  gas  at  the 
output  of  the  gas  separator  contains  nitrogen,  carbon  dioxide, 
water  vapor,  and  methanol  vapor.  After  reducing  the  water  content 
by  the  condenser  at  room  temperature,  the  gas  temperature,  the 
pressure,  and  the  carbon  dioxide  volume  concentration  of  the 
water  saturated  gas  can  be  measured. 

At  side  B  air  is  used  as  gas  which  can  be  humidified  before 
entering  the  test  cell.  Because  of  the  illustrated  processes  inside  the 
cell  (see  Fig.  2X),  the  output  of  the  cathode  consists  of  a  mixture  of 
air,  water,  and  carbon  dioxide.  The  content  of  these  gases  can  be 
determined  by  the  humidity  sensor,  the  carbon  dioxide  sensor,  and 
the  pressure  sensor. 

The  DMFC  operating  mode  is  used  to  investigate  the  correlations 
between  cell  voltage,  methanol  permeation,  and  water  output  at 
the  cathode  as  a  function  of  operating  parameters  such  as  current 
density,  methanol  concentration,  volume  flow  at  anode  and  cath¬ 
ode,  and  cell  temperature. 

2.2.2.  Pseudo-half-cell  operating  mode 

The  phc  mode  is  very  similar  to  the  DMFC  mode.  The  only  dif¬ 
ference  is  that  the  gas  at  the  cathode  side  is  nitrogen  and  therefore 
the  cell  needs  electrical  energy  to  oxidize  the  methanol  at  the 
anode.  In  this  mode  only  two  processes  inside  the  MEA  are  related 
to  the  production  of  carbon  dioxide:  the  methanol  oxidation  at  the 
anode,  see  process  1  in  Fig.  2Y,  which  generates  the  carbon  dioxide 
in  the  anode  flow  field,  and  the  carbon  dioxide  diffusion  from 
anode  to  cathode,  see  process  2  in  Fig.  2Y.  The  strength  of  the 
carbon  dioxide  diffusion  from  anode  to  cathode  can  be  determined 
by  measuring  the  carbon  dioxide  content  in  the  output  of  the 
cathode. 


X 

DMFC  mode 

side  A  |  side  B 


side  A  |  side  B 


input  side  B: 

N2,  h2o 


Fig.  2.  Illustration  of  the  different  processes  inside  the  MEA  in  DMFC  mode  (X)  and  phc  mode  (Y). 
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This  mode  is  used  to  get  an  approximate  idea  about  the  anode 
polarization  curve  as  a  function  of  the  operating  parameters.  The 
permeated  methanol  cannot  be  oxidized  at  the  cathode  so  that  the 
measurement  of  the  carbon  dioxide  concentration  at  the  cathode 
output  provides  a  maximum  value  for  the  diffusion  of  carbon  di¬ 
oxide  from  anode  to  cathode  during  DMFC  operation.  In  DMFC 
mode  the  diffusion  of  carbon  dioxide  should  be  smaller  because  the 
gradient  in  the  carbon  dioxide  concentration  between  anode  and 
cathode  is  smaller  when  carbon  dioxide  is  produced  at  the  cathode 
side,  too. 

2.2.3.  Reverse  operating  mode 

Apart  from  the  polarity  of  the  power  source,  the  reverse  oper¬ 
ating  mode  is  identical  to  the  phc  mode.  The  reverse  mode  uses  the 
methanol  permeation  through  the  MEA  where  methanol  is 
oxidized  at  side  B  (anode)  forming  carbon  dioxide,  which  can  be 
measured  by  the  carbon  dioxide  sensor  B.  At  the  cathode  (side  A), 
hydrogen  is  produced,  see  process  1  in  Fig.  3X.  The  produced  carbon 
dioxide  can  diffuse  from  anode  to  cathode,  see  process  2  in  Fig.  3X. 

The  reverse  mode  offers  the  possibility  for  a  voltage  based 
approximation  of  the  methanol  permeation.  By  increasing  the 
voltage  the  current  density  will  be  limited  because  the  permeation 
of  the  methanol  to  the  cathode  is  limited  by  the  physical  properties 
of  the  MEA.  The  current  density  of  the  methanol  permeation  jp  is 
approximated  by  the  measured  limiting  current  density  jum  [!]• 

2.2.4.  Calibration  operating  mode 

The  calibration  method  presented  in  this  study  is  a  modification 
of  the  reverse  operating  mode.  The  correlation  between  the  electric 
current  and  the  produced  quantity  of  carbon  dioxide  depends  on 
the  quantity  of  side  products  of  the  methanol  oxidation. 

In  case  that  an  unknown  part  of  the  methanol  is  oxidized  to 
formic  acid  or  methanal,  the  correlation  between  the  electric  cur¬ 
rent  and  the  quantity  of  carbon  dioxide  is  not  clear.  In  order  to 
generate  carbon  dioxide  with  a  known  correlation  between  the 
electric  current  and  the  quantity  of  carbon  dioxide,  formic  acid  is 
electrochemically  oxidized  instead  of  methanol  because  this  reac¬ 
tion  has  been  investigated  intensively  and  is  free  from  side  re¬ 
actions  [8-10].  Equation  (1)  shows  the  simplified  electrochemical 
oxidation  process  for  the  formic  acid. 

CHOOH  C02  +  2H+  +  2e“  ( 1 ) 


The  formic  acid  permeates  through  the  MEA  and  can  be  oxidized 
on  side  B  of  the  MEA,  see  process  1  in  Fig.  3Y.  The  carbon  dioxide  is 
mixed  with  the  nitrogen  within  the  flow  field  and  a  mixture  of 
carbon  dioxide,  nitrogen,  and  water  flows  out  of  the  cell.  If  it  is 
assumed  that  the  gas  is  saturated  with  water  and  the  gas  temper¬ 
ature  and  the  ambient  pressure  are  measured,  the  gas  composition 
can  be  calculated.  To  avoid  that  carbon  dioxide  diffuses  from  side  B 
to  side  A,  the  formic  acid  solution  must  be  saturated  with  carbon 
dioxide.  For  this  reason,  the  valve  at  side  A  in  Fig.  1  directs  the  liquid 
phase  back  to  the  tank  so  that  the  solution  is  pumped  in  a  circle  and 
is  saturated  after  a  time  with  carbon  dioxide.  To  minimize  the  time 
for  saturation,  a  small  volume  of  100  ml  high  concentrated  aqueous 
formic  acid  solution  in  combination  with  a  high  flow  rate  of 
12  ml  min-1  for  the  formic  acid  solution  is  used.  The  volume  of  the 
flow  field  channel  is  1  ml.  Therefore,  the  retention  time  of  the  so¬ 
lution  in  the  test  cell  is  approximately  5  s. 

A  further  reason  for  a  high  concentrated  formic  acid  solution  is 
the  low  permeation  rate  in  comparison  to  methanol.  According  to 
Rhee,  the  permeation  of  formic  acid  through  a  Nafion®  membrane 
reaches  its  maximum  at  concentrations  about  10  mol  l-1  [11], 
which  is  why  this  concentration  is  used  for  the  calibration  exper¬ 
iments  in  this  study. 

To  avoid  a  decomposition  of  formic  acid  at  the  platinum  elec¬ 
trodes  to  hydrogen  and  carbon  dioxide  at  higher  temperatures,  the 
calibration  is  carried  out  at  room  temperature.  Hence,  all  pre¬ 
heaters  in  the  experimental  setup,  see  Fig.  1,  are  removed.  At  side  B, 
the  moisture  sensor  is  not  required  for  the  calibration  as  well  as  the 
condenser.  For  the  calibration  process  the  carbon  dioxide  sensor  is 
located  behind  the  pressure  sensor  at  side  B  in  the  experimental 
setup.  By  means  of  the  temperature  sensor  in  front  of  the  pressure 
sensor,  the  temperature  of  the  output  gas  can  be  controlled. 
Although  the  electrochemical  decomposition  of  the  formic  acid 
produces  heat,  no  increase  in  gas  temperature  is  observed.  Obvi¬ 
ously,  the  high  flow  rate  of  the  formic  acid  solution  at  side  A  in 
combination  with  the  cooling  effect  by  water  evaporation  at  side  B 
compensate  the  heating. 

2.3.  Calculation  of  the  gas  composition 

The  composition  of  the  gases  at  the  cell  outputs  of  side  A  and 
side  B  in  Fig.  1  depends  on  the  operating  mode.  In  this  section,  the 
gas  composition  for  the  calibration  mode  and  the  related 


side  A  |  side  B 


input  side  B: 

N2,  h2o 


output  side  A: 
formic  acid,  <— 
H20,  C02 

Formic  acid 
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i 

power 


A  a  A  a  A 

flow  fields  — j — • — |— 1 

electrodes . j . . . ; 

membrane - - 


input  side  B: 
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Fig.  3.  Illustration  of  the  different  processes  inside  the  MEA  in  reverse  mode  (X)  and  calibration  mode  (Y). 
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uncertainty  is  calculated.  The  gas  compositions  for  the  DMFC  and 
phc  mode  are  determined  in  Appendix  A.  The  calculations  assume 
that  all  gases  behave  like  an  ideal  gas  and  that  the  gas  which  passes 
the  condenser  is  saturated  with  water  at  temperature  Tx  (X  =  A  for 
side  A,  X  =  B  for  side  B).  Therefore,  the  temperature  Tx  is  the  dew 
point  of  the  gases. 


2.3.1.  Calibration  mode 

The  theoretical  value  of  the  carbon  dioxide  volume 
concentration  a  is  a  function  of  the  experimental  quantities  such  as 
the  electric  current  J,  the  nitrogen  flow  rate  frn,  the  gas  temperature 
T,  and  the  ambient  pressure  pa. 

According  to  Equation  (1)  and  Faraday’s  law,  the  amount  of 
carbon  dioxide  per  time  n'Cd  can  be  described  as  a  function  of  the 
electric  current  I.  Relating  to  standard  conditions  the  carbon  di¬ 
oxide  flow  rate  frCd  can  be  calculated  by  means  of  the  general  gas 
constant  R  (8.314  J  mol-1  I<-1),  the  Faraday  constant  F 
(96,485  C  mol-1),  the  standard  temperature  Ts  (273.15  I<),  and  the 
standard  pressure  ps  (101325  Pa). 

n'cd  =  60  /  (2  Fr1  mol  min-1  (2) 

frcd(/)  =  60  /  (2  FT1  ■RTs(psy'lm3  min-1  1 

For  an  ideal  gas,  the  ratio  of  the  water  flow  rate  (frw)  and  the 
flow  rates  of  the  dry  gases  (frn  -  nitrogen,  frCd  -  carbon  dioxide)  is 
equal  to  the  ratio  of  the  partial  pressures.  Under  the  condition  that 
the  total  pressure  is  equal  to  the  ambient  pressure  pa,  this  yields: 

frw- (frn  +  frcd)  1  —  Pw'(Pn+Pcd)  1  —  Pw'(Pa— Pw)  1  (3) 

pw  -  partial  pressure  of  water,  pn  -  partial  pressure  of  nitrogen,  pCd 
-  partial  pressure  of  carbon  dioxide. 

The  partial  pressure  of  water  is  a  function  of  gas  temperature  T. 
All  measurements  presented  in  this  study  were  performed  at  room 
temperature  (23  °C)  and  ambient  pressure  pa.  Furthermore  the  gas 
is  assumed  to  be  saturated  with  water.  To  calculate  the  saturated 
water  vapor  pressure  pw  on  the  basis  of  the  temperature  T,  the 
parameters  of  the  Magnus  formula  12].  A,  m,  and  Tc  were  fitted  to 
the  partial  pressure  data  of  water  in  the  range  between  10  °C  and 
40  °C.  The  data  were  downloaded  from  Ref.  [13]  and  the  software 
originPro  was  used  to  fit  the  parameters. 

Pw  =  Pw(T)  =  AM-exp(m-T-{T  +  Tcy^ 

Fit  results  :  AM  =  (610.68  +  /  -  0.03)  Pa  (4) 

m  =  (17.307  +  / -0.004) 

Tc  =  (237.72  +  / -0.07)  °C 

Although  the  value  for  the  accuracy  of  the  temperature  trans¬ 
mitter  is  0.4  °C,  it  has  to  be  considered  that  the  temperature  T  in 
Equation  (4)  corresponds  to  the  dew  point  temperature  of  the  gas. 
To  take  a  small  difference  between  the  gas  temperature  and  the 
dew  point  temperature  into  account,  the  accuracy  for  the  deter¬ 
mination  of  the  dew  point  temperature  is  estimated  at  ±2  °C. 
Because  of  the  low  relative  accuracy  of  T  in  comparison  to  the 
values  for  the  fit  parameters  in  Equation  (4)  the  function  pw  can  be 
assumed  to  be  only  a  function  of  the  random  variable  T.  According 
to  the  Gaussian  error  propagation,  the  standard  deviation  of  pw  is: 

dpw  =  m  Tc  pw  (Tc  +  r)“2  dT;dT  =  2  °C  (5) 

In  line  with  Equations  (2)  and  (3),  the  carbon  dioxide  volume 
concentration  ob  in  the  gas  mixture  behind  the  condenser  at  side  B 
in  Fig.  1  can  be  calculated  as  follows  (unless  otherwise  noted  for 
calibration  is  a3  =  a): 


V  =  frcd  •  (frn  +  frcd  +  frw)  1 

frw  =Pw'(Pb— Pw)  '  (frn  +  frcd ) 

->•  O  =  frcd  •  (PB  -  Pw)  ■  ((frn  +  frcd)  'Pa)  1 
-  ff  =  /•(pB-Pw)'(Pa'(/  +  C'frn))-1;c  =  2  F  Ps-(60  R  Tsy1 

(6) 

a  is  a  function  of  the  random  variables  /,  frn,  pa,  and  pw.  These 
variables  are  assumed  to  be  statistically  independent.  Therefore, 
the  standard  deviation  of  a  can  be  determined  by  Gaussian  error 
propagation.  The  standard  deviations  of  these  values  are  estimated 
by  the  accuracy  of  the  metering  devices  used.  Based  on  this 
approach,  the  standard  deviation  of  a  can  be  calculated  by  Equation 
(7),  where  dx  describes  the  standard  deviation  of  the  measured 
unit  x. 


da  =  (tf  +  t|  +  t|  +  t|)“°'5 

t\  =  da- (5/)_1  dl  da-(dl)~^ 

t2  =  a(7-(afrn)_1-dfrn  a<7-(afrn)“ 

t3  =  a(7-(apa)_1  -dpa  a<7-(apa)“ 

£4  —  a<j"(apw)  '  dpw  a<j'(apw) 


-1 

1 

-1 


a  ’C-  frn  •(/•(/  +  c*  frn))  1 
=  ^(J  +  C-frnr1 

-  &’Pw’ (Pa’ (Pa’Pw)) 

—  ~<J'  ( Pa  'Pw) 


(7) 


Equations  (6)  and  (7)  enable  the  carbon  dioxide  volume  con¬ 
centration  a  and  its  uncertainty  to  be  calculated  for  each  adjust¬ 
ment  of  I  and  frn  if  the  gas  temperature  T  and  the  ambient  pressure 
pa  are  known. 


2.4.  Reference  gases 

In  order  to  verify  the  calibration  quality,  commercial  calibration 
gases  (Linde)  with  different  carbon  dioxide  contents  are  used. 
Table  1  shows  the  gas  concentrations  employed  for  carbon  dioxide 
in  nitrogen. 


2.5.  Basic  inspections 

The  manufacturer  of  the  sensor  states  the  insensibility  to  water 
vapor  [14],  only  the  dilution  of  the  gas  by  water  vapor  has  to  be 
considered.  To  prove  this  statement,  the  calibrated  sensor  is  flushed 
with  a  commercial  calibration  gas.  The  measured  value  is  compared 
with  the  value  for  humidified  gas.  If  the  gas  is  saturated  with  water 
vapor  and  the  gas  temperature  is  known,  the  dilution  effect  can  be 
calculated  and  compared  with  the  difference  of  the  measured 
values.  The  calibration  gas  with  the  highest  concentration  of  19.9% 
v/v  is  used  for  this  test,  because  with  this  gas  the  sensor  has  the 
highest  relative  accuracy,  see  Section  3.5. 

For  the  purpose  of  demonstrating  the  performance  of  the  cali¬ 
bration  method,  a  carbon  dioxide  sensor  with  a  large  measuring 
range  of  20%  v/v  is  used.  According  to  the  manufacturer,  the 


Table  1 

Carbon  dioxide  concentrations  employed  for  verifying  the  calibration.  The  uncer¬ 
tainty  is  related  to  a  confidence  level  of  95%. 


C02  volume  concentration  in  %  v/v 

Uncertainty  in  %  v/v 

0.203 

0.002 

0.999 

0.010 

1.515 

0.015 

2.00 

0.02 

3.00 

0.03 

3.570 

0.036 

5.00 

0.05 

10.0 

0.1 

14.900 

0.149 

19.900 

0.199 
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accuracy  of  this  sensor  at  25  °C  is  ±(1.5%  of  metering  range  ±2%  of 
reading),  including  repeatability,  nonlinearity,  and  calibration  un¬ 
certainty.  This  applies  for  concentrations  above  2%  of  full  scale.  As  a 
result,  for  a  metering  range  of  20%  v/v  the  signal  must  be  higher 
than  0.4%  v/v.  At  this  point,  it  has  to  be  noted  that  the  carbon  di¬ 
oxide  measurement  device  consists  of  two  parts:  the  transmitter 
(GMT220)  and  the  sensor  (GMP221 ).  Connecting  and  disconnecting 
the  sensor  from  the  transmitter  can  change  the  reading  value  so  if 
the  transmitter  is  changed,  the  calibration  has  to  be  repeated.  In 
this  paper,  the  sensor  in  combination  with  a  transmitter  is  cali¬ 
brated  as  one  device.  This  is  one  of  the  reasons  for  the  improvement 
in  accuracy. 

The  signal  noise  and  the  repeatability  of  the  signals  have  to  be 
determined  to  identify  the  accuracy  attainable  with  the  sensor 
using  the  method  presented.  In  order  to  quantify  the  level  of  signal 
noise  and  signal  stability  within  the  metering  range,  the  sensor  is 
flushed  with  nitrogen  for  several  days  and  the  signal  is  recorded. 
Additionally,  the  noise  level  of  the  sensor  for  a  carbon  dioxide 
concentration  of  20%  v/v  is  determined  by  flushing  the  sensor  with 
commercial  carbon  dioxide  calibration  gas.  The  measured  standard 
deviation  of  the  signals  can  be  used  to  estimate  the  signal  noise  as  a 
function  of  the  reading. 

In  a  further  test,  the  repeatability  of  the  signals  has  to  be 
checked.  The  sensor  is  calibrated  initially  and  afterward  the  cali¬ 
bration  is  verified  by  measuring  the  carbon  dioxide  concentration 
of  the  commercial  calibration  gases.  The  sensor  is  then  turned  off. 
Two  weeks  later,  the  verification  is  repeated  in  order  to  quantify  the 
repeatability  of  the  results. 

3.  Results  and  discussion 

In  this  section,  the  results  of  the  preliminary  tests  are  presented. 
An  error  analysis  is  used  to  determine  the  necessary  adjustments 
for  the  electric  current  I  and  the  nitrogen  flow  rate  frn.  Afterward, 
the  calibration  procedure  is  explained  and  the  measured  data  are 
presented.  Polynomial  functions  are  fitted  to  the  data  and  are  used 
to  correct  the  measured  raw  data  for  the  calibration.  The  accuracy 
of  the  calibration  is  verified  by  the  investigation  of  commercial 
calibration  gases  with  a  known  carbon  dioxide  concentration. 
Finally,  the  statement  of  the  manufacturer  about  the  insensibility  of 
the  carbon  dioxide  sensor  to  water  vapor  is  proven. 

3.1  Sensor  noise 


a/%  v/v 

Fig.  4.  Approximation  of  the  standard  deviation  as  a  function  of  the  carbon  dioxide 
concentration.  The  measured  data  are  approximated  by  a  linear  sensor  noise  function 
sn.  The  parameters  of  sn  are  calculated  by  the  noise  level  of  a  =  0%  v/v  and  a  =  20%  v/v. 

3.2.  Error  analysis 

All  random  variables  considered  in  this  study  are  assumed  to  be 
normally  distributed.  Gaussian  error  propagation  is  used  to  calcu¬ 
late  the  complete  uncertainty  of  Equations  (7)  and  (8)  because  of 
the  statistical  independence  of  the  data  from  the  carbon  dioxide 
transmitter,  the  electric  load,  and  the  mass  flow  meter  for  the  ni¬ 
trogen  flow  rate. 

do  total  =  ((da)2  +  sn2)  °'5  (9) 

To  calibrate  the  sensor,  different  adjustments  for  the  electric 
current  I  and  the  flow  rate  frn  are  calculated  by  means  of  Equation 
(6)  to  generate  a  carbon  dioxide-nitrogen  mixture  with  carbon 
dioxide  volume  concentrations  a  between  0%  v/v  and  20%  v/v.  The 
metering  range  of  the  mass  flow  meter  has  a  maximum  value  for 
the  flow  rate  of  100  ml  min-1.  In  addition,  the  formic  acid  perme¬ 
ation  through  the  MEA  is  restricted.  The  absolute  maximum  current 
in  this  study  is  0.9  A.  The  calculated  values  are  presented  in  Fig.  5. 

In  order  to  illustrate  the  different  contributions  to  the  calculated 
accuracy,  Fig.  6  shows  the  single  values  of  each  term  of  the 


The  data  of  the  basic  inspections  (see  Section  2.5)  are  used  to 
estimate  the  signal  noise  of  the  sensor  by  calculating  the  standard 
deviation  of  the  measured  data.  This  value  may  be  different  for  each 
sensor.  A  value  of  5  10-4  was  determined  for  the  employed  device 
using  nitrogen  as  the  purge  gas.  To  estimate  the  noise  level  as  a 
function  of  the  carbon  dioxide  concentration,  the  sensor  is  flushed 
with  commercial  calibration  gases  for  30  min.  The  concentrations 
used  are  listed  in  Table  1.  The  standard  deviations  of  the  measured 
data,  presented  in  Fig.  4,  show  an  almost  linear  dependency  of  the 
noise  level  as  a  function  of  a.  To  demonstrate  that  a  low  concen¬ 
tration  (nitrogen)  and  a  high  concentration  (20%  v/v)  measurement 
are  sufficient  to  estimate  this  function,  the  function  parameters  are 
calculated  by  using  these  two  noise  values. 

sn(cr)  =  0.00256  ±0.00228  a  (8) 

The  resulting  two-point  linearization  is  shown  in  Fig.  4. 

In  summary,  it  follows  that  the  calibration  method  presented  in 
this  paper  needs  only  one  calibration  gas  with  a  high  concentration 
to  estimate  the  noise  level  function.  This  function  is  needed  for  the 
error  analysis  in  the  following  section. 


a  /  %  v/v 

Fig.  5.  Presentation  of  the  different  adjustments  for  the  electric  current  I  and  the  ni¬ 
trogen  flow  rate  frn  for  the  generation  of  a  carbon  dioxide  mixture  with  a  carbon  di¬ 
oxide  volume  concentration  a. 
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Fig.  6.  Contributions  to  the  standard  deviation  of  <7totai-  A  fitted  second  degree  poly¬ 
nomial  function  shows  the  total  standard  deviation  as  a  function  of  a. 


calculated  standard  deviation  including  the  signal  noise,  see  also 
Equation  (9). 

According  to  Fig.  6,  the  uncertainty  of  the  mass  flow  meter  is 
responsible  for  the  largest  contribution  to  the  standard  deviation, 
followed  by  the  uncertainty  of  the  ampere  meter.  The  uncertainty 
of  the  water  partial  pressure  provides  the  smallest  contribution  and 
the  uncertainty  of  the  ambient  pressure  determination  is  negli¬ 
gible.  The  standard  deviation  is  dominated  by  the  signal  noise  sn  of 
the  sensor  for  values  smaller  than  1%  v/v.  A  second  degree  poly¬ 
nomial  function  is  fitted  to  the  total  values  of  the  standard  devia¬ 
tion  to  describe  the  total  standard  deviation  as  a  function  of  a. 


function  is  used  to  estimate  the  uncertainty  of  the  calibrated  values 
ac  (see  Equations  (11)  and  (12)). 

3.3.  Calibration  procedure 

To  avoid  a  systematic  error,  the  order  of  the  adjustments  pre¬ 
sented  in  Fig.  5  is  permuted  for  the  calibration  procedure.  After 
each  setting  of  frn  and  I  measurements  are  only  made  after  the 
measured  signal  is  stable.  Additionally,  the  measured  pressure  pa 
and  gas  temperature  T  is  sent  to  the  transmitter  to  update  the  in¬ 
ternal  signal  correction.  Then  the  sensor  values  are  recorded  for 
10  min.  The  mean  and  standard  deviation  of  the  measured  data  are 
calculated  for  each  adjustment. 

To  demonstrate  the  nonlinear  behavior  of  the  sensor,  the 
theoretical  carbon  dioxide  values,  according  to  Equation  (6),  are 
plotted  in  Fig.  7  A  as  a  function  of  the  calculated  mean  values.  A 
linear  function  is  fitted  to  the  data  and  the  residuals  are  plotted  in 
Fig.  7B.  This  residual  plot  shows  that  a  part  of  the  sensor  signal 
cannot  be  described  by  a  linear  function.  Additionally,  the  dotted 
lines  in  Fig.  7B  show  that  the  expected  uncertainty,  according  to 
Equation  (10),  is  smaller  than  the  residuals. 

For  a  further  reduction  of  the  residuals,  a  third  degree  poly¬ 
nomial  is  chosen  to  fit  the  data  for  values  smaller  than  or  equal  to 
1.5%  and  a  fifth  degree  polynomial  is  used  to  fit  the  data  for  values 
higher  than  1.5%.  The  software  originPro  is  used  for  the  fit  process. 
The  following  equations  show  the  calibration  function  ac  as  a 
function  of  the  measured  value  of  a. 

ac  =  a0  +  ar(7  +  a2-<72  +  a3-<73  *  0%  v/v  <  a  <  \.S%  v/v 
a0  =  -0.00148;  aq  =  0.569;  a2  =  27.3;  a3  =  -539 

(11) 


(jc  =  a0  +  ai  -  o'  +  a2-o-2  +  a3  -  a3  +  a4-(j4  +  a5-(j5  ;  1.5  %  v/v  <  a  <  20%  v/v 

Qq  =  —0.00404;  Gq  =  1.05;  a2  —  —1.02;  d3  -  —13.6;  U4  =  145;  Q5  =  —319 


da total  fit  =  0.0035  +  0.0073  (7  +  0.0002  (72 


(10) 


This  function  quantifies  the  limits  of  the  calibration  technique 
due  to  the  devices  employed.  In  the  following  sections,  this 


The  residuals  for  the  fitted  polynomial  functions  (presented  in 
Fig.  8)  are  smaller  than  the  expected  uncertainties  according  to 
Equation  (10). 


0 


•  calculated  values  as  a  function  of  the  measured  data 
- fitted  linear  function 


2^ 


£  0.0 


•  residuals 

- positive  values  according  to  Equation  10 

. negative  values  according  to  Equation  10 


10  15 

a  /  %  v/v 


10  15 

o7  %  v/v 


Fig.  7.  A)  a  represents  the  mean  values  of  the  calibration  procedure  and  the  calculated  values  of  the  theoretical  data,  according  to  Equation  (6).  A  linear  function  is  fitted  to  the  data. 
B)  The  differences  between  the  theoretical  data  and  the  linear  function  are  shown  in  the  residual  plot.  The  dotted  lines  show  the  expected  uncertainty  of  the  calibration. 
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a  I  %  v/v  a  /  %  v/v 


Fig.  8.  A)  Residuals  for  Equation  (11)  and  concentrations  smaller  than  1.5%  v/v.  B)  Residuals  for  Equation  (12)  and  concentrations  higher  than  1.5%  v/v.  All  residuals  are  smaller  than 
the  expected  uncertainty,  according  to  Equation  (9). 


At  the  end  of  this  section,  it  can  be  concluded  that  the  calibra¬ 
tion  functions  (Equations  (11)  and  (12))  can  be  used  to  reduce  the 
residuals  of  the  fit  functions  so  that  the  expected  uncertainties, 
according  to  Equation  (10),  are  higher  than  the  residuals.  This 
represents  the  limit  of  the  presented  calibration  method  because 
this  uncertainty  is  determined  by  the  sensor  noise  and  the  devices 
used  for  this  calibration.  The  quality  of  this  calibration  is  verified  in 
the  next  section. 

3.4.  Verification  of  the  calibration 

For  the  following  verification  procedure  the  calibration  func¬ 
tions  Equations  (11 )  and  (12)  are  used  to  correct  the  sensor  signals. 
Equation  (10)  is  applied  to  estimate  the  uncertainties.  The  sensor  is 
flushed  with  each  of  the  commercial  calibration  gases  (see  Table  1 ) 
until  no  signal  change  within  the  scattering  is  measurable.  After¬ 
ward,  the  corrected  signal  is  recorded  for  10  min.  The  data  are  used 
to  calculate  the  mean  and  the  standard  deviation  (Equation  (10))  as 
an  estimator  of  the  true  value  and  the  uncertainty  of  the  mea¬ 
surement.  To  show  the  results  more  clearly,  the  differences  be¬ 
tween  the  means  and  the  given  values  for  the  carbon  dioxide 
concentrations  are  considered.  In  this  presentation,  the  set-point 
values  are  zero  with  the  uncertainties  specified  in  Table  1.  Fig.  9 
shows  the  measured  values  (mean(<7c)  -  ^specified)  where  the  er¬ 
ror  bars  indicate  the  uncertainty,  which  is  estimated  by  twice  the 
standard  deviation  (Equation  (10)).  Based  on  the  assumption  of 
normally  distributed  and  statistically  independent  variables,  the 
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Fig.  9.  Experimental  verification  of  the  calibration.  The  unfilled  diamonds  show  the 
difference  between  the  measured  values  and  the  specified  values  of  the  commercial 
calibration  gas.  All  specified  values  are  within  the  error  bars  of  the  measured  values. 


error  bars  in  Fig.  9  approximate  the  range  to  the  confidence  level  of 
95%.  The  specified  values  indicate  that  these  error  bars  are  reliable 
estimators  of  the  uncertainty  of  the  measurements  because  all 
values  presented  as  black  points  in  Fig.  9  are  within  these  error 
bars. 

This  verification  experiment  was  repeated  after  two  weeks.  No 
differences  could  be  measured  within  the  uncertainties. 

The  relative  accuracy  of  the  sensor  can  be  determined  by  means 
of  Equation  (10).  Fig.  10  shows  the  calculated  relative  accuracy  as  a 
function  of  the  corrected  value  for  the  carbon  dioxide  concentra¬ 
tion  in  comparison  to  the  manufacturer’s  specifications.  In  both 
cases  a  confidence  level  of  95%  is  chosen. 

Fig.  10  clearly  shows  the  improvement  of  sensor  accuracy  by  the 
calibration.  For  concentrations  higher  than  0.2%  v/v  the  accuracy 
was  enhanced  up  to  a  factor  of  30  in  comparison  with  the  manu¬ 
facturer’s  specifications.  Within  the  range  of  1%  v/v  and  20%  v/v  the 
relative  accuracy  of  the  calibrated  sensor  is  2%  (confidence 
level  =  95%). 

3.5.  Influence  of  water  vapor  on  the  sensor  signal 

To  prove  the  influence  of  water  vapor  on  the  signal  of  the  carbon 
dioxide  sensor,  the  sensor  is  flushed  with  commercial  calibration 
gas  with  19.9%  v/v  carbon  dioxide  and  the  signal  is  recorded.  Then, 


ac  /  %  v/v 


Fig.  10.  Relative  accuracy  of  a  carbon  dioxide  measurement  (confidence  level  =  95%). 
Metering  device:  GMT220  Series  Carbon  Dioxide  Transmitters  (Vaisala),  metering 
range  of  the  sensor  (Type:  GMP221 ):  20%  v/v.  Calculated  relative  accuracy  of  the  sensor 
reading  (see  Equation  (10))  after  calibration  in  comparison  with  the  manufacturer’s 
specifications. 
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the  gas  is  humidified  with  water  at  room  temperature  and  the 
temperature  of  the  gas  is  measured.  Assuming  that  the  gas  is 
saturated  with  water,  the  saturated  water  vapor  pressure  can  be 
determined  by  Equation  (4).  The  carbon  dioxide  volume  concen¬ 
tration  ahum  of  the  humidified  gas  can  be  calculated  by  the  defi¬ 
nition  of  a  in  combination  with  Equation  (3). 


confirmed,  only  the  dilution  effect  of  water  on  the  gas  concentra¬ 
tion  has  to  be  considered  for  the  measurements.  In  this  study,  the 
dilution  effect  of  water  is  considered  in  all  calculations  of  gas 
compositions. 


ffhum  =  frCd  •  (frn  +  frCd  +  frw)  1 ;  Eq.  3  :  —  frw  =  Pw(TB)  •  (pa  -  Pw)  1  ■  (frn  +  frCd) 
°hum  =  frcd-((frn  +  frCd)-(l +Pw(7'b)-(Pb-Pw(7'b))“1^) 


(13) 


The  parameters  Ohum,  Pw(7b),  Pb»  and  a dry  are  statistically  inde¬ 
pendent,  therefore  the  uncertainty  can  be  calculated  by  Gaussian 
error  propagation. 


dA 

dA 


theor  =  |^dry)  +  (^hum 
theor  =  A  theor  •  ^  (dp  w  ( TB )  ’  (Pw  (Tb)) 

d^dry  ’  (°dry) 


*(Pb)  j)2  (14) 


According  to  Equation  (13),  the  difference  between  the  carbon 
dioxide  volume  concentration  of  the  dry  gas  and  the  humidified  gas 
(Atheor)  can  be  calculated  by  the  concentration  of  the  dry  gas 
multiplied  by  the  factor  PwC^bMPb)-1  which  can  be  determined  by 
the  gas  temperature  TB  in  combination  with  Equation  (4)  and  the 
ambient  pressure  pB.  This  value  can  be  compared  with  the  exper¬ 
imental  determined  concentration  difference  Aexp  =  0b, dry  -  0b, hum- 
The  results  of  the  measurement  are  shown  in  Fig.  11. 

The  temperature  of  the  humidified  reference  gas  is  23  °C.  Ac¬ 
cording  to  Equation  (4),  the  water  vapor  pressure  pw(TB)  is  2811  Pa. 
The  ambient  pressure  is  pa  =  98,280  Pa.  The  result  of  Equation  (13) 
and  Equation  (14)  for  Atheor  is  (0.57  ±  0.03)%  v/v.  Within  the  un¬ 
certainty  of  the  sensor  values  this  value  cannot  be  distinguished 
from  the  measured  value  of  (0.38  ±  0.34)%  v/v.  Therefore,  the 
statement  of  the  carbon  dioxide  sensor  manufacturer  can  be 


Influence  of  water  vapor  on  the  carbon  dioxide  sensor  signal 
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Fig.  11.  Measurement  of  the  carbon  dioxide  volume  concentration  with  and  without 
humidification.  The  gas  temperature  and  the  ambient  pressure  were  measured  to 
calculate  the  theoretical  difference  between  the  humidified  and  dry  gas 
concentrations. 


4.  Relevance  of  the  calibration  for  direct  methanol  fuel  cells 
investigations 

In  this  section,  it  is  demonstrated  how  the  research  results  of 
direct  methanol  fuel  cell  efficiency  can  be  improved  by  means  of 
the  calibrated  carbon  dioxide  sensor.  For  this  purpose,  the  cali¬ 
brated  sensor  is  used  to  detect  the  signal  0a,  see  experimental  setup 
in  Fig.  1.  A  second  carbon  dioxide  sensor  (same  type,  metering 
range:  5%  v/v)  is  calibrated  with  the  presented  method  and  used  in 
the  setup  to  measure  0-B.  The  MEA  which  is  employed  for  the 
following  experiments  is  of  the  same  type  as  the  MEA  used  for 
calibration  (see  Section  2.1). 

4.1.  Objective  target  of  the  DMFC  investigations 

One  objective  target  in  investigating  direct  methanol  fuel  cell 
efficiency  is  to  increase  the  Faraday  efficiency  ?7f,  which  is  defined 
by  Equation  (15). 

%  =j-(j+jpf1  (15) 

A  typical  approach  to  increase  pp  is  to  reduce  the  methanol 
permeation  by  using  new  materials  for  the  MEA  preparation. 
Another  option  is  to  choose  appropriate  operation  conditions,  so 
that  the  permeation  current  density  is  as  small  as  possible.  The 
latter  requires  an  understanding  about  the  correlation  between  jp 
and  the  operating  parameters  such  as  electric  current  density  j, 
methanol  concentration  c(MeOH),  volume  flow  at  the  cathode  vfB, 
cell  temperature,  and  anodic  volume  flow  vfA.  During  the  variation 
of  the  operational  parameters  the  permeation  current  density  jP 
and  the  diffusion  current  density  jo  will  change.  The  permeation 
current  density  can  be  determined  on  the  base  of  the  carbon  di¬ 
oxide  signal  <tb  in  case  the  diffusion  current  density  is  known,  see 
Equation  (35). 

4.2.  Determining  the  diffusion  current  density  jo  in  DMFC  mode 

The  diffusion  current  density  can  be  calculated  by  means  of  the 
carbon  dioxide  signal  a  a  (see  Equation  (36))  if  the  nitrogen  injector 
of  the  gas  separator  (see  Fig.  1,  side  A)  removes  as  much  carbon 
dioxide  from  the  liquid  phase  of  the  anode  outlet  that  the  dissolved 
carbon  dioxide  in  the  liquid  phase  is  negligible  within  the  uncer¬ 
tainty  of  the  measurement. 

The  principle  is  that  the  higher  the  volume  flow  the  better  the 
nitrogen  removes  the  carbon  dioxide  from  the  liquid  phase  but  the 
smaller  is  the  carbon  dioxide  volume  concentration  (jB  and  thus  the 
accuracy  of  the  signal.  For  this  reason,  the  volume  flow  must  be  as 
small  as  possible  but  as  high  as  necessary.  To  find  the  best  adjust¬ 
ment  for  the  nitrogen  injector,  the  experimental  setup  is  driven  in 


156 


M.  Stabler,  A.  Burdzik  /  Journal  of  Power  Sources  262  (2014)  147-161 


B 


E 

o 

< 


0.025 

0.024 

0.023 

0.022 

0.021 

0.020 

0.019 

0.018 

0.017 

0.016 

0.015 

0.014 

0.013 

0.012 

0.011 

0.010 


time  /  s 


time  /  s 


CQ 

< 

t> 


Fig.  12.  A):  Variation  of  the  injected  nitrogen  volume  flow  vfn  (unit:  ml  (min  cm2)-1)  into  the  gas  separator  (see  Fig.  1,  side  A)  during  constant  operation  conditions.  Operational 
parameters:)  =  0.1  A  cm-2,  c(MeOH)  =  0.75  mol  l-1,  anodic  volume  flow  0.57  ml  (min  cm2)-1,  cathodic  volume  flow  11  ml  (min  cm2)-1,  and  cell  temperature  of  70  °C.  B):  Same  data 
as  in  A  but  the  uncertainties  are  calculated  by  using  the  manufacturer’s  specification  (da  =  0.5  (0.015  metering  range  +  0.02  reading  value),  metering  range:  5%  v/v).  Error  bars  are 
calculated  by  Equations  (37)  and  (43)  and  show  single  standard  deviation.  Additionally,  the  measured  carbon  dioxide  concentration  of  aA  and  (7B  are  plotted. 


phc  mode.  For  a  non-zero  electric  current  density  j ,  the  operating 
parameters  of  the  MEA  are  fixed  and  the  volume  flow  of  the  ni¬ 
trogen  injector  vfn  at  the  gas  separator  of  the  anode  side  is  varied. 
For  each  adjustment  of  vfn  the  carbon  dioxide  signal  aA  is  used  to 
calculate  j’a  (see  Equation  (29))  and  the  difference  j  -  jA  which  is 
compared  with  the  diffusion  current  density  jo.  Jd  is  calculated  by 
means  of  the  measured  carbon  dioxide  concentration  <7B,  the 
pressure  pe,  and  the  gas  temperature  Tb  at  the  cathode,  see  Equa¬ 
tion  (42).  If  the  nitrogen  removes  all  carbon  dioxide  from  the  outlet 
of  the  anode,  a  further  increase  of  the  nitrogen  flow  will  not  change 
the  value  of  j  -  j’a.  The  results  for  such  a  measurement  are  shown  in 
Fig.  12. 

The  nitrogen  volume  flow  vfn  is  varied  between  5.7  and 
28.3  ml  (min  cm2)-1  during  constant  operation  conditions  for  the 
MEA.  Fig.  12A  shows  that  the  increase  from  22.7  to 
28.3  ml  (min  cm2)'1  does  not  result  in  a  further  change  of  j  -  j’a.  The 
values  for  j  -jA  and  jo  are  comparable  within  the  uncertainty  of  the 
measurement.  The  best  adjustment  of  the  nitrogen  volume  flow  for 
the  chosen  operating  parameters  is  22.7  ml  (min  cm2)-1.  Fig.  12B 
shows  the  same  data  but  the  uncertainties  of  the  values  are 
calculated  by  using  the  manufacturer’s  specifications  for  the 


Fig.  13.  Required  nitrogen  volume  flow  to  remove  the  carbon  dioxide  from  the  liquid 
phase  of  the  anodic  outlet  as  a  function  of  the  electric  current  density. 


uncertainty  of  the  carbon  dioxide  sensor.  Additionally,  the 
measured  carbon  dioxide  concentrations  are  plotted  in  Fig.  12B  to 
show  that  the  values  differ  between  a  few  percent  and  0.25%  v/v. 
The  improved  calibration  of  the  carbon  dioxide  sensor  over  a  wide 
metering  range  allows  a  more  precisely  adjustment  of  the  nitrogen 
volume  flow. 

The  appropriate  nitrogen  volume  flow  is  determined  with  the 
presented  method  for  different  electric  current  densities.  The 
anodic  volume  flow  is  constantly  kept  at  0.57  ml  (min  cm2)-1  for  all 
measurements.  The  determined  values  for  vfn  as  a  function  of  the 
current  density  are  presented  in  Fig.  13.  A  power  function  is  fitted  to 
the  data  to  calculate  the  required  nitrogen  flow  rate  on  the  base  of 
the  adjusted  electric  current  density  j. 

The  function  in  Fig.  13  is  used  in  the  following  experiments  to 
adjust  the  nitrogen  volume  flow  for  the  gas  separator  as  a  function 
of  the  electric  current  density  j.  To  prove  this  relation,  a  voltage 
current  curve  is  measured  in  phc  mode  with  several  adjustments 
for  j  and  a  methanol  concentration  of  0.5  mol  1_1.  The  comparison  of 
the  data  in  Fig.  14  shows  that  the  measured  diffusion  current 
density  jo  is  within  the  uncertainty  of  the  calculated  diffusion 
current  density  jD  =  j  -jA. 

With  the  method  presented  in  this  section  the  diffusion  current 
density  jo  can  be  determined  by  the  difference  of  j  and  jA  in  phc 
mode.  This  method  can  also  be  used  in  DMFC  mode  because  the 
operation  of  the  anode  side  is  identical  in  both  modes.  The  option  to 


j  /A- cm-2 

Fig.  14.  Result  from  the  measured  voltage  current  curve  in  phc  mode.  Difference)  -)A 
and  )D  as  a  function  of  the  electric  current. 
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j/A-cnr2  j/A-cnr2 

Fig.  15.  Measurement  of  the  methanol  permeation  and  the  diffusion  current  density  as  a  function  of  the  electric  current  density  for  two  different  methanol  concentrations. 


calculate  the  diffusion  current  density  jo  in  DMFC  mode  enables  the 
calculation  of  the  permeation  current  density  jP,  according  to 
Equation  (35). 

4.3.  Optimizing  the  Faraday  efficiency 

One  strategy  for  the  optimization  of  the  Faraday  efficiency  is  to 
analyze  the  correlation  between  the  permeation  current  density  jP 
and  the  operating  parameters  to  find  adjustments  with  a  high  ef¬ 
ficiency.  As  an  example,  Fig.  15  shows  the  measured  values  for  jP 
and  jo  as  a  function  of  the  electric  current  density  j  and  for  two 
different  methanol  concentrations.  Fig.  15A  shows  the  values 
measured  in  the  DMFC  mode.  The  nitrogen  volume  flow  for  the  gas 
separator  on  side  A  (see  Fig.  1 )  is  adjusted,  according  to  Section  4.2, 
as  a  function  of  the  electric  current  density.  Equations  (36)  and  (37) 
are  used  to  calculated  jo  =j  -  Ja  and  the  corresponding  uncertainty. 
The  values  for  jP  are  calculated  by  Equations  (35)  and  (38)  is 
employed  to  determine  the  uncertainty.  The  error  bars  in  Fig.  15 A 
show  the  single  standard  deviations  based  on  the  presented  cali¬ 
bration  method  and  Fig.  15B  shows  the  same  data  with  single 
standard  deviations,  calculated  by  means  of  the  manufacturer’s 
specification  (da  =  0.5  (0.015  metering  range  +  0.02  reading), 
metering  range:  5%  v/v). 

According  to  the  data  presented  in  Fig.  15,  the  uncertainty  for 
the  determination  of  the  permeation  and  diffusion  current  density 
is  reduced  by  the  employed  calibration  technique  about  a  factor  of 
three  to  four.  This  improvement  factor  depends  on  the  operating 
parameters  used  for  this  measurement.  For  example,  the  smaller 
the  volume  flow  vfp  at  the  cathode  of  the  MEA  the  smaller  the 
uncertainty  of  jP  in  Fig.  15A  (see  Equation  (38)  influence  of  term  ti). 

To  quantify  the  effect  of  this  improvement  on  the  Faraday  effi¬ 
ciency,  the  error  propagation  of  p?  as  a  function  of  j  and  jP  has  to  be 
determined.  Fig.  15  shows  that  jP  correlates  with  j,  therefore  the 
Gaussian  error  propagation  cannot  be  used.  According  to  Brandt 
[15],  the  error  propagation  of  correlated  units  has  to  be  calculated 
by  transforming  the  covariance  matrix  of  the  units  by  means  of  the 
Jacobian  matrix  J(f)  of  the  analyzed  function/.  Equation  (16)  shows 
the  Jacobian  matrix  of  p^. 

J(Vf )  =  [a??F-(9/r\  3^F-(9/p)_1] 

=  [jp  •  0'  +  jp)  2 ,  -J •  0'  +  jp)  2]  (16) 

The  transformation  of  the  covariance  matrix  [var (j),  cov(jjP); 
cov(jjP),  var(jP)]  by  the  Jacobian  matrix  of  pF  provides: 


(cfyp)2  =  (j+jp)  4-((jp  djf  +  (j-djpf  -2  •j-jp-covO'jp)) 

(17) 

To  find  an  estimation  for  co v(jjp)jp  is  approximated  by  a  linear 
function  of  j.  Because  of  the  linear  behavior  of  the  covariance 
function,  this  provides: 

jP  =  a  +  bj 

->cov(j,jP)  =  cov(j,  a  +  b  -j)  =  co  v(j,a)  (18) 

+  b-cov(j.j)  =  bvar(j)  =  bdj2 

The  term  co v(j,a)  is  equal  to  zero  because  there  is  no  correlation 
between  the  current  density  j  and  the  constant  a.  Equation  (19) 
quantifies  the  uncertainty  of  p?  as  a  function  of  j,  jP,  dj ,  djP,  and  b. 

(drip)2  =  (j+jPy4-{(jp-dj)2  +  (j-djp)2 

-2  j  jp  b  dj2^  ;for  djp  see  Equation  (35)  (19) 

The  factor  b  is  calculated  by  fitting  a  linear  function  to  the 
permeation  current  density.  For  a  methanol  concentration  of 
0.75  mol  1_1  b  is  equal  to  -0.163  and  for  0.50  mol  l-1  b  is  equal 
to  -0.174.  With  Equation  (15),  Equation  (19),  and  the  values  for  b 
the  Faraday  efficiency  for  the  data  in  Fig.  15  can  be  calculated. 

Fig.  16  shows  the  Faraday  efficiency  p?  as  a  function  of  j  for  two 
methanol  concentrations.  In  Fig.  16A  the  uncertainty  dp^, a  is 
calculated,  according  to  Equation  (7),  in  combination  with  Equation 
(35)  and  da  of  the  calibrated  sensors.  In  Fig.  16B  the  manufacturer’s 
specification  is  used  to  calculate  dp p,b  (Equation  (7)  in  combination 
with  Equation  (35)  and  da  =  0.5  (0.015  metering 
range  +  0.02  reading),  metering  range:  5%  v/v).  The  ratio 
d?7F,B-(d?7F,A)-1  is  plotted  in  Fig.  16B.  The  uncertainty  of  the  Faraday 
efficiency  is  reduced  by  a  factor  of  about  3  to  4  by  the  use  of  the 
presented  calibration  method. 

Additionally,  Fig.  16  shows  the  efficiency  pf  calculated  for  the 
case  that  all  carbon  dioxide  in  the  cathode  outlet  is  assumed  to  be 
formed  by  methanol  permeation.  If  the  carbon  dioxide  is  not  being 
considered,  the  calculated  efficiency  is  smaller  than  the  real 
Faraday  efficiency.  This  fact  can  be  more  clearly  distinguished  by 
the  measured  data  in  case  that  the  presented  calibration  method  is 
used.  In  case  the  manufacturer’s  specification  is  applied  to  calculate 
the  uncertainty  of  p' F,  the  difference  between  p F  and  pf  is  within 
the  single  standard  deviation. 

Summarizing  the  above,  it  can  be  said,  that  the  presented  cali¬ 
bration  method  is  able  to  reduce  the  uncertainty  of  the  determined 
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Fig.  16.  Faraday  efficiency  ^  as  a  function  of  electric  current  density  for  two  different  methanol  concentrations.  Further  operating  parameters:  temperature  70  °C,  anode  volume 
flow  0.57  ml  (min  cm2)-1,  cathode  volume  flow  (20  +  50-j)  ml  (min  cm2)-1  for  j  <  100  mA  cm-2;  (10  +  250-j)  ml  (min  cm2)-1  for  j  >  100  mA  cm-2.  A):  Data  with  calculated 
uncertainties,  according  to  Equation  (7).  B):  Same  efficiency  data  as  A  with  calculated  uncertainties,  according  to  Equation  (7),  in  combination  with  the  manufacturer’s  specification 
for  the  uncertainties  of  the  measured  carbon  dioxide  concentration.  Additionally,  the  lines  show  the  ratio  of  the  uncertainties  of  B  and  A,  calculated  for  each  point  of  j. 


permeation  current  density  jp  and  the  diffusion  current  density  jD 
about  a  factor  of  up  to  4.  The  uncertainty  of  the  Faraday  efficiency 
can  be  reduced  by  a  factor  of  4. 


5.  Conclusions 

A  calibration  technique  for  carbon  dioxide  sensors  is  introduced 
in  the  presented  paper.  The  basic  idea  of  this  procedure  is  to  employ 
the  available  devices  of  a  DMFC  test  rig  in  order  to  generate  a 
mixture  of  carbon  dioxide  and  nitrogen  as  calibration  gas.  Carbon 
dioxide  sensors  are  often  used  in  DMFC  test  stations  to  analyze  the 
methanol  permeation  in  MEAs.  If  a  relative  precision  of  a  few 
percent  is  sufficient,  they  can  be  used  as  an  alternative  to  expensive 
and  complex  analyzing  methods  such  as  FTIR  or  mass  spectroscopy. 
Carbon  dioxide  sensors  can  be  easily  integrated  in  automated  test 
rigs  and  they  are  well  suited  for  real-time  measurements.  A 
disadvantage  is  the  low  accuracy  which  depends  strongly  on  the 
sensor’s  measurement  range.  Furthermore,  they  have  to  be  cali¬ 
brated  regularly. 

In  this  study,  the  necessary  preliminary  tests  and  error  analysis 
are  presented  in  detail.  The  calibration  method  and  the  subsequent 
verification  of  the  calibration  by  the  use  of  commercial  calibration 
gases  are  described.  It  is  shown  that  the  relative  accuracy  of  the 
calibrated  sensor  can  be  improved  by  a  factor  of  up  to  30, 
depending  on  the  metering  range.  Additionally,  the  relative  sensor 
accuracy  is  constant  (2%,  confidence  level  =  95%)  within  95%  of  the 
full  scale. 

It  is  demonstrated  how  the  improvement  of  the  calibration 
quality  enhanced  the  accuracy  of  the  determined  Faraday  efficiency 
in  direct  methanol  fuel  cell  investigations.  By  means  of  the  cali¬ 
brated  sensors,  it  is  shown  that  the  diffusion  current  density  in 
DMFC  operating  mode  can  be  identified  more  precisely  in  com¬ 
parison  to  the  situation  that  the  manufacturer’s  specification  about 
the  accuracy  is  used  to  find  out  the  uncertainty  of  the  measure¬ 
ments.  Because  of  a  smaller  uncertainty  of  the  Faraday  efficiency  by 
using  the  presented  calibration  technique,  the  influence  of  the 
carbon  dioxide  diffusion  from  anode  to  cathode  on  the  determi¬ 
nation  of  ?7f  is  clearly  visible. 

The  demonstrated  technique  is  suitable  for  a  precise,  fast,  and 
fully  automatic  calibration  of  carbon  dioxide  sensors. 


Appendix  A 

This  appendix  will  provide  the  formulary  to  calculate  the  gas 
composition  of  the  gases  behind  the  condensers  in  the  experi¬ 
mental  setup,  illustrated  in  Fig.  1,  for  the  two  operating  modes: 

1.  DMFC  mode 

2.  phc  mode 

Additionally,  the  uncertainty  for  each  formula  is  determined. 

In  comparison  to  the  calculations  of  the  gas  composition  in 
calibration  mode,  see  Section  2.3,  the  measurements  of  the  modes 
described  here  provide  the  carbon  dioxide  concentrations.  In  these 
cases  the  operator  is  interested  in  the  permeation  current  density  jp 
and  the  diffusion  current  density  jo  which  are  defined  in  this 
section. 

As  a  reminder:  The  calculations  assume  that  all  gases  behave 
like  an  ideal  gas  and  that  the  gas  which  passed  the  condenser  is 
saturated  with  water  at  temperature  Tx  (X  =  A  for  side  A,  X  =  B  for 
side  B,  see  Fig.  1). 

1 )  Gas  composition  for  DMFC  operating  mode 

Side  A 

According  to  Fig.  1,  when  the  gas  has  passed  the  condenser  the 
amount  of  substance  n\  depends  on  the  flow  rate  of  nitrogen  frn 
injected  into  the  gas  separator,  the  flow  rate  of  carbon  dioxide 
frj(C02)  which  is  correlated  to  the  adjusted  electric  current  density 
j,  the  flow  rate  of  carbon  dioxide  hpC CO2)  which  diffuses  from 
anode  to  cathode,  expressed  as  diffusion  current  density  jD,  and  the 
flow  rate  of  water  frw. 


n'A  =  frn  +  frj(C02)  -  frjD(C02)  +  frw  (20) 

The  operating  parameter  for  the  mass  flow  controller  is  the  flow 
rate.  All  flow  rates  in  this  study  are  related  to  standard  conditions: 
ps  =  101325  Pa,  Ts  =  273.15  K.  Therefore,  the  flow  rate  of  nitrogen 
can  be  described  as  a  function  of  volume  flow  vf  A,  the  general  gas 
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constant  R  (8.314  J  mol  1  K  1),  the  standard  temperature  Ts,  and 
standard  pressure  ps. 

frn  =  (Ji-Ts)-1-ps-vf,A  (21) 

The  flow  rate  of  carbon  dioxide  fr,(C02),  correlating  with  the 
adjusted  electric  current  density  j,  can  be  described,  according  to 
Equation  (1)  and  Faraday’s  law,  as: 

frj(C02)  =  j-A-(6-F)-'l-,A  =  17.6  cm2  is  the  MEA  area. 

(22) 

The  flow  rate  of  carbon  dioxide  which  diffuses  from  anode  to 
cathode  can  be  expressed  as  a  diffusion  current  density  jD: 


djl  =  (ti  )2  +  (f2)2  +  (t3)2 

ti  =  ajA-(0<7A)“1-d<rA  =  (1  -/A)-jA-d( ta-((1  -/a -<ta)-<ta)“1 
t2  =  0JA  •  (0V;)-1  •  dvfA  =  jA  •  dv  f  A  •  (vfA)_1 
h  =  9jA  (0/A)_1d/A  #=  jA  dvfA  (l  -/a  +  <7a)_1 
dfl  =  (dPw(PA)"1)  +  (dPA'Pw(PA)"2) 

(30) 

The  parameter  dpw  is  calculated,  according  to  Equation  (5). 

Equations  (29)  and  (30)  allow  the  description  of  the  measured 
carbon  dioxide  at  side  A  as  a  current  density  jA  which  can  be 
compared  with  the  adjusted  electrical  current  density  j  to  validate 
the  setup  and  the  calculations. 


frjD(C02)  =  jD-A-(6-F)_1  (23) 

It  is  not  possible  to  differentiate  frj(C02)  and  AptCCb)  by  the 
measurement  of  one  carbon  dioxide  signal,  therefore  these  two 
flow  rates  are  combined  to  frCd- 

frcd  =  frj(C02)  -  frjD(C02)  =  <j-jD)-A-{6-F)~' 

=  jA-A-(6-F)~'-,jA=j-jD  (24) 

For  ideal  gases,  the  ratio  of  water  flow  rate  frw  to  the  total  dry 
gas  flow  rate  can  be  expressed  as  a  ratio  of  the  partial  pressure  of 
water  to  the  partial  pressure  of  gases  without  water,  see  Equation 
(3).  Therefore,  the  flow  rate  of  water  can  be  calculated  by  Equation 

(25). 

frw  —  Pw'(Pa—  Pw)  1  '  (frn  +  frcd)  (25) 

By  means  of  Equations  (21)— (25).  Equation  (20)  can  be  written 
in  the  following  form: 

riA  =  ((R-Is)1  'Ps'Vf;  +jA-A-(6-F)-1)  -(pA-(pA  —  p(H20)))_1 

(26) 

Equation  (27)  shows  the  definition  of  the  measured  carbon  di¬ 
oxide  volume  concentration. 

*a  =  frcd -Kr1  (27) 

Insertion  of  Equations  (23)  and  (26)  provides: 

aA  =  (jA -A  -  (6-F)-1 )  •  ((R-  Ps)”1  -ps  •  VA 

+ Ja -A- (6-F)"1)  •  (pA-  (pA  -  p(H20))_1)  1  (28) 

After  a  few  transposition  steps  the  formula  for  jA  is  obtained. 

Ja  =  <rA-g-vfA-(l  -/a-^a)-1;  vfA  =  vf^-A'1 

g  =  PsF-6(RTs)~' 

/a  =/(Ta)  =Pw(Pa)-1  =  •  exp  (m  •  TA  •  (TA  +  Tc)_1 )  •  (pA)-1 

(29) 

The  variables  in  Equation  (29)  of  jA  are  statistically  independent 
because  crA,  vfA,  and  the  variables  for  /A  are  measured  with  inde¬ 
pendent  devices.  The  uncertainty  djA  is  calculated  by  Gaussian  error 
propagation. 


Side  B 

The  calculation  of  the  gas  composition  at  the  output  of  the 
condenser  of  side  B  is  carried  out  analog  to  the  section  above.  The 
gas  flow  consists  of  air  (fra),  injected  at  the  input  of  side  B.  A  part  of 
the  oxygen  is  reduced  by  the  electrochemical  reaction  (fr,(02))  and 
the  oxidation  of  the  permeated  methanol  (frJP(02)).  Additionally, 
the  gas  flow  contains  carbon  dioxide  (frJP(C02))  from  the  oxidation 
of  the  permeated  methanol  and  the  diffusion  from  anode  to  cath¬ 
ode  (frjo(C02)).  This  defines  the  equation  for  the  gas  flow  rate  at  the 
output  of  the  condenser  of  side  B  (n'B). 

nB  =  fra  -  frj  (02)  -  fijp  (02) + frjp  (C02 )  +  frjD  (C02) + frw  (31) 

The  flow  rate  of  air  can  be  described  as  a  function  of  volume  flow 
vf  b,  analog  to  Equation  (21).  With  the  following  relations  Equation 
(31)  can  be  transformed. 

fra  =  (J?-rs)_1-ps-vf'B 
frj(02)  =  j-A-(4-F)-1 
frjp(02)  =JP-A-(4-F)-1 

frjP(C02)  =jp-A-(6-F)-1  (32) 

frjo(C02)  =jD-A-(6-F)-1 

frw  =  Pw(7b) •  (PB  —  Pw(Jb))_1  •  (fra  -  frj(02)  -  frjp(02) 

+  frjp(C02)  +  frjD(C02)) 

After  a  few  transposition  steps  Equation  (33)  is  obtained. 
n'B  =  (i?Ts)-1-ps-vf;+AF-1-(jp-6-1  +JD-6-1  -Jp-4-1 
-Jd'4_1)  •Pw(7'b)'(Pb  -Pw(Tb))-1 

(33) 

Analog  to  Equation  (27),  the  definition  of  the  carbon  dioxide 
volume  concentration  at  the  output  of  the  condenser  at  side  B  (o-B) 
is  given  by  Equation  (34). 

(7B  =  frcdK)”1  (34) 

By  inserting  Equation  (33)  into  Equation  (34)  and  isolating  jp  the 
methanol  permeation  current  density  can  be  described  as: 

Jp  =  (or- (g  vfB  -  6-j-4_1) 

_i  (35) 

-  JV  (1  -  ^b -/(Tb)))  •  (l  +  ^b -2-1 -/(Tb)^ 

For  the  definition  of  the  constant  g  and  the  function  /  see 
Equation  (29);  vfB  =  vfB-A_1. 


160 


M  Stabler,  A.  Burdzik  /  Journal  of  Power  Sources  262  (2014)  147-161 


For  the  calculation  of  the  permeation  current  density  jp  the 
diffusion  current  density  j'd  is  required.  By  means  of  Equation  (29), 
the  current  density  j'a  can  be  calculated  on  the  base  of  the  measured 
carbon  dioxide  concentration  a  a.  Assuming  that  all  carbon  dioxide 
can  be  removed  from  the  solution  in  the  gas  separator  (see  side  A  in 
Fig.  1 )  the  diffusion  current  density  results  from  the  difference  of 
the  electric  current  density  j  and  j'a. 

J'd  =  j  -Ja  (36) 

Because  j  and  j’a  are  statistically  independent,  the  uncertainty  of 
j'd  can  be  determined  by  Gaussian  error  propagation. 

djl  =  ( dj )2  +  (djA)2  (37) 

With  Equation  (37)  the  uncertainty  of  jp  can  be  calculated  by 
Gaussian  error  propagation,  too. 

4/p  =  (fl)2  +  (f2)2  +  (f3)2  +  (f4)2  +  (f5)2 

c=l+0.5-<7B-/(TB) 

ti  =  djp  •  (a<jB) 

*=  (1  -/(Tb))- fe-vfB-6-4-1  -O+Jd)) -d<rB-c-2 

t2  =  0Jp  •  (SVb)-1  •  dvfB  =  g- •  dvfB  •  c-1 

f3  =  0jp '  ( dj)  dj  =  6  •  4  1  (Tgdj-c^1  (38) 

f4  =  9/p  •  (9Jd)_1  '  djo  —  (1  -  <^b  — /  (7b))  •  djD  •  c-1 
ts  =  djP-(df(Tjl)rt-df(TB) 

=  (6  •  4-1  •  ffB  -jD + ffB  •  (g  •  vf  B  -  6  •  4- 1  •  j)  •  c-2 )  •  d/(7B) 
d/(TB)2  =  (dpw •  (p(7'b))“1  ) 2  +  (dp(TB)  Pw •  (P(7b))“2)2 

Equations  (35)  and  (36)  allow  to  calculate  the  permeation  cur¬ 
rent  density  jp  and  the  diffusion  current  density  jo  for  the  DMFC 
operation  mode. 

2)  Gas  composition  for  phc  operating  mode 


Side  A 

The  gas  composition  at  side  A  in  phc  mode  is  equal  to  the  sit¬ 
uation  at  side  A  in  DMFC  mode. 


Side  B 

The  gas  flow  injected  at  the  input  of  side  B  consists  of  nitrogen 
(frn).  Additionally,  the  gas  flow  contains  hydrogen  (frj(H2))  produced 
by  the  electrochemical  reaction  and  carbon  dioxide  which  diffuses 
from  anode  to  cathode  (frjo(C02)).  This  defines  the  equation  for  the 
gas  flow  rate  at  the  output  of  the  condenser  of  side  B  (n'e). 

n'B  =  frn  +  frj(H2)  +  frjD(C02)  +  frw  (39) 

The  flow  rate  of  air  can  be  described  as  a  function  of  volume  flow 
vf  b,  analog  to  Equation  (21).  With  the  following  relations  Equation 
(39)  can  be  transformed. 

fra  =  (RTs)-:psvf'B 
frj(H2)  =  j-A-(2-F)~J 
frjD(C02)  *  jD-A-(6-F)  1 

frw  =  Pw(Tb)-(Pb  —  Pw(TB))_1  •  (fra  +  frj(H2)  +  frjD(C02)) 

(40) 

Insertion  of  these  expressions  in  Equation  (39)  provides: 


n'„  =  ((P-rs)“1ps-vfB 

+  A  ■  F1  •  (j ■  2-1  +  jo  ■  6-1 )  )  -Pw(Tb)  •  (PB  -  Pw(rB))"1 

(41) 

The  definition  of  the  carbon  dioxide  volume  concentration  at 
the  output  of  the  condenser  at  side  B  (ctb)  is  given  by  Equation  (34). 
Insertion  of  Equation  (41 )  yields  an  expression  to  determine  jo. 

Jd  =  (ffB-(g-vfB+6-j-2-1))-(l-ffB-/(TB)r1  (42) 

For  the  definition  of  the  constant  g  and  the  function  /  see 
Equation  (29);  vf  B  =  vfB  A  1 

The  uncertainty  of  jo  can  be  calculated  by: 

dj2o  =  (fi)2  +  (f2)2  +  (f3)2  +  (t4)2 

k  =  1  -  ffB  ~f(TB) 

ti  =  ajD-(a(TBr1-d(TB  =  (1  -f(TB))-  (g’vfs  +  6-2-1  -j)  •  d<rB - k— 2 
h  =  ajo  ■  fai'/g)  1- dvfs  =  g- (JB  -c/i/b  ■  k  1 
t3  =  ajD-(0j)_1  dj  =  6-4 <rBdjk-i 

t4  =  ajo^a/^u))-1 -d/(TB)  =  <TB(g-vfB  +  6-2-^j)-k~2-df(TB) 

df(TB)2  =  (dpw-(p(TB))-^2  +  (dp(TB)-pw-(p(TB)r2)2 

(43) 

The  diffusion  current  density  jo  in  Equation  (42)  provides  an 
opportunity  to  prove  the  experimental  setup.  The  measurement  of 
j'a  enables  the  calculation  of  the  diffusion  current  density  jD  by 
means  of  the  difference  j  -  j'a  (Equation  (36))  which  can  be 
compared  with  the  measured  carbon  dioxide  concentration  in 
Equation  (42).  These  values  must  be  the  same. 

In  summary  of  Appendix  A,  the  presented  equations  allow  the 
calculation  of  the  permeation  current  density  jp  and  the  diffusion 
current  density  j'd  on  the  base  of  the  operational  parameters  j,  vfA, 
vfB,  (7a,  (7 B,  Pa,  Pb,  Ta,  and  TB.  The  accuracy  of  the  carbon  dioxide 
sensor  influences  the  expression  t\  in  Equations  (30)  and  (38). 
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List  of  symbols 

A:  active  area  of  the  MEA 

frw:  water  flow  rate 

frcci:  carbon  dioxide  flow  rate 


M.  Stabler,  A.  Burdzik  /  Journal  of  Power  Sources  262  (2014)  147-161 


161 


frn:  nitrogen  flow  rate 

vfA.B-  volume  flow  at  side  A  or  B 

v/a.B-'  volume  flow  at  side  A  or  B  standardized  to  MEA  area 
F:  Faraday  constant  (96,485  C  mol-1) 

I:  electric  current 
T:  gas  temperature 

a:  carbon  dioxide  volume  concentration 

ac:  corrected  carbon  dioxide  volume  concentration 

a specified •'  carbon  dioxide  concentration  of  commercial  calibration  gas 

dutotai •'  total  standard  deviation  of  the  carbon  dioxide  volume  concentration 

pCd:  partial  pressure  of  carbon  dioxide 

pw:  saturated  water  vapor  pressure 


pn:  partial  pressure  of  nitrogen 

pa:  ambient  pressure 

ps:  standard  pressure  (101,325  Pa) 

Ts:  standard  temperature  (273.15  K) 

R:  general  gas  constant  (8.314  J  K-1  mol-1) 

Am,  rn,  Tc:  Magnus  parameters 

sn:  signal  noise  function 

j:  electric  current  density 

jp:  permeation  current  density 

jD:  carbon  dioxide  diffusion  current  density 

g:  constant:  g  =  ps-F-6-(l?Ts)_1 


